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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC 
(SI)  UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

inches 

feet 

square  feet 

cubic  feet  per  second 

feet  per  second 
per  second 

degrees  (angular) 

Fahrenheit  degrees 


$2 

25.4 

0.3048 

0.09290304 

0.02831685 

0.3048 

0.01745329 

5/9 


To  Obtain 

millimetres 

metres 

square  metres 

cubic  metres  per  second 

metres  per  second 
per  second 

radians 

Celsius  degrees  or 
Kelvins* 


* To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F) 
readings,  use  the  following  formula:  C ■ (5/9) (F  - 32).  To  obtain 

Kelvin  (K)  readings,  use:  K ■ (5/9) (F  - 32)  - 273.15. 
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CONVEX  BOUNDARIES  IN  CONVERGING  SUPERCRITICAL  FLOW 


PART  I:  INTRODUCTION 

Statement  of  the  Problem 


1.  Steep  chutes  are  commonly  used  in  hydraulic  structures  as 
conveyances  for  high-velocity  flows  between  a spillway  crest  and  an 
energy  dissipator.  Normally,  the  type,  width,  and  elevation  of  the 
crest  are  determined  by  the  stage-discharge  relationship  prescribed  by 
overall  project  requirements.  Similarly,  the  width  and  elevation  of  the 
floor  at  the  entrance  to  the  energy  dissipator  (often  a stilling  basin) 
are  determined  by  the  need  to  match  the  hydraulic  characteristics  of  the 
dissipator  with  the  stage-discharge  relationship  prescribed  by  down- 
stream channel  requirements.  Whenever  the  crest  length  is  greater  than 
the  width  of  the  energy  dissipator,  the  chute  must  converge  in  the  down- 
stream direction.  The  subject  of  this  study  is  the  shape  in  plan  of  the 
chute  sidewalls;  the  object  is  to  study  the  relationships  commonly  used 
to  design  an  acceptable  rate  of  curvature  of  the  sidewalls.  The  rate  of 
curvature  is  a function  of  the  hydraulic  characteristics  of  the  flow  as 
well  as  of  the  geometrical  characteristics  imposed  on  the  chute  by 
general  project  requirements. 

Practical  Application 

2.  The  project  features  that  most  commonly  enter  into  the  overall 
chute  spillway  design  are  the  approach  channel,  wing  walls,  spillway 
crest,  chute,  and  energy  dissipator  (Figures  la  and  lb).  In  hydraulic 
engineering  manuals  and  in  research  studies,  these  items  are  often 
treated  separately;  in  the  overall  project,  however,  the  designs  must 

be  Integrated  so  that  acceptable  flow  conditions  exist  throughout. 
Hydraulic  model  tests  are  the  most  common  design  tool  used  to  select  the 
best  overall  geometry  in  situations  not  previously  encountered  and  not 
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STEEP  CONVERGING 


ELEVATION  or  ENERGY  GRADE  LINE  OF 
APPROACH  FLO* 

ELEVATION  OF  WATER  SURFACE  OF 
APPROACH  FLO* 

= VELOCITY  HEAD  OF  THE  APPROACH  FLO* 
= CREST  HEIGHT 

SLOPE  OF  THE  CHUTE 

incremental  distance  islopei 

INCREMENTAL  DISTANCE  (HORIZONTAL! 

= INCREMENTAL  WIDTH  (HORIZONTAL! 


• PLAN 


b PROFILE 


Figure  1.  Typical  steep  chute 

amenable  to  analytical  solution.  A few  general  considerations  that  in- 
fluenced the  approach  taken  in  this  study  are  outlined  in  the  following 
paragraphs  for  each  of  the  above  project  features. 

3.  The  upstream  approach  to  the  spillway  should  align  the  flow 
normal  to  the  crest — an  improperly  designed  approach  can  impose  flow 
disturbances  that  are  amplified  along  the  chute  and  are  unacceptable  in 
the  downstream  channel. The  wing  walls  usually  should  extend  upstream 
from  the  crest  so  that  the  initial  contraction  occurs  in  a subcrltlcal 
flow  regime;  they  may  be  curved  in  plan  so  as  to  prevent  vortex  shedding 
at  the  upstream  edge. 
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4.  The  design  procedure  for  the  profile  of  the  spillway  crest  is 
well  covered  in  hydraulic  engineering  technical  literature;  in  particu- 
lar, design  charts  are  available  in  Corps  of  Engineers  Hydraulic  Design 

A1 

Criteria  (HDC).  The  crest  geometry,  in  plan,  may  be  straight  or  curved 

as  indicated  by  the  solid  and  broken  lines,  respectively,  in  Figure  la. 

The  crest  may  be  gated;  but  more  commonly,  no  gates  are  present  and  the 

crest  alone  controls  the  flow  (overfall  spillway).  Flow  disturbances 

originating  at  piers  and  abutments  on  the  crest  may  cause  substantial 

D2 

surface  irregularities  in  the  downstream  channel. 

5.  One  of  the  most  common  causes  of  spillway  failures  has  been 
the  improper  design  of  steep  chutes;  criteria  for  the  crest  shape  are 
well  known  but  the  requirement  for  the  correct  shaping  of  the  chute 
sidewalls  in  plan  is  often  overlooked.  Straight  walls,  shown  as  solid 
lines  in  Figure  1,  are  economical  and  perform  satisfactorily  in  some 
physical  situations;  however,  these  walls  may  cause  unacceptable  surface 
waves  in  the  downstream  channel.  The  remedy  is  to  use  gently  curving 
sidewalls  as  indicated  previously  by  the  dashed  lines  in  Figure  la.  The 
curvature  of  the  sidewalls  then  becomes  a delicate  problem  since  the 
flow  in  the  chute  is  supercritical,  and  whenever  this  high-velocity  flow 
encounters  an  improperly  aligned  sidewall,  large  dynamic  phenomena  are 
created.  Obvious  consequences  are  standing  waves,  overtopping,  and  ex- 
cessive structural  loadings  as  well  as  a nonuniform  flow  in  the  down- 
stream channel. 

6.  Energy  dissipation  below  the  chute  is  normally  accomplished  by 
means  of  plunge  pools,  deflector  buckets,  or  stilling  basins.  For  plunge 
pools  and  deflector  buckets  some  degree  of  nonuniformity  may  be  accept- 
able in  the  flow  as  it  leaves  the  chute,  depending  on  the  physical  char- 
acteristics of  the  exit  channel.  However,  in  the  case  of  a stilling 
basin,  the  flow  should  be  as  nearly  uniform  as  possible. 
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Review  of  the  Literature 


7.  Converging  chutes  have  been  used  in  numerous  projects  and 
often  successfully.  The  high  percentage  of  successes  is  due  to  the 
major  role  that  physical  hydraulic  models  play  in  a comprehensive  hy- 
draulics design  program.  Grand  Coulee  Dam,^  for  example,  which  uses  a 
parallel-walled  chute,  required  14  separate  models  representing  designs 
from  the  planning  to  the  construction  phase  in  order  to  obtain  the  final 
adopted  configuration.  Surface  shock  waves  in  the  flow  in  the  chute  are 
easily  observed  in  a reasonably  sized  spillway  model.  Consequently,  the 
configuration  of  the  chute  (or  appurtenant  structural  components)  can  be 
modified  in  the  model  to  eliminate  (or  reduce  to  acceptable  levels)  the 
amplitude  of  the  waves.  Of  course,  in  order  to  reduce  the  effort  in  the 
laboratory,  a good  initial  design  is  a tremendous  advantage. 

8.  Common  configurations  of  chute  spillways  and  the  associated 

design  criteria,  as  of  31  March  1965,  of  the  Corps  of  Engineers  are 

A2 

given  in  Engineer  Manual  (EM)  1110-2-1603  ; comparable  information  for 

flood  control  channels  as  of  1 July  1970  is  contained  in  EM  1110-2- 
A3 

1601.  The  present  suggested  criteria,  including  data  from  nine  exist- 

C1 

ing  chutes,  were  reviewed  by  R.  G.  Cox  and  published  as  Engineer 

A4 

Technical  Letter  (ETL)  1110-2-158  during  1972.  Design  criteria  have 

A5 

also  been  described  for  the  Bureau  of  Reclamation  and  in  numerous  hy- 
draulic engineering  treatises  (References  A6,  A7,  and  A8,  for  example). 
An  overview  of  published  design  guidance  for  situations  in  which  stand- 
ing waves  are  to  be  avoided  and  uniform  flow  is  required  at  the  entrance 

to  the  downstream  energy  dissipator  is  as  follows. 

A4 

a.  Optimum  flow  conditions  occur  with  a crest  formed  by  a 
break  in  invert  grade  or  by  a low  sill  formed  as  an  inte- 
gral part  of  the  chute  slope;  this  precludes  ogee  crests 
with  toe  curves,  although  this  type  of  design  has  been 
used  successfully  in  practice.  Curving  the  chute  crest 
does  not  appreciably  affect  the  flow  conditions  in  the 
chute . 

b.  Straight,  parallel  sidewalls  are  preferred  whereas  reverse 
curves  are  normally  unacceptable;  where  required, 
gradually  convergent  or  divergent  walls  can  be  used. 
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A4 

c.  Vertical  walls  are  preferred  ; if  sloping  walls  are  used, 
the  side  slope  should  be  considered  relative  to  the  chute 
invert  rather  than  to  the  horizontal  plane. 

ci.  The  walls  should  be  straight  in  the  vicinity  of  the  crest; 
the  convergence  ratio  at  the  crest,  |As/Az|,  should  be 
greater  than  or  equal  to  5.0.  The  incremental  distance 
along  the  chute,  As,  and  the  incremental  width,  Az,  are 
shown  in  Figure  1.  The  straight  walls  should  extend  down- 
stream beyond  the  point  where  the  local  Froude  number  be- 
comes greater  than  1.5.A4  The  local  Froude  number. 


/gd 


is  based  on  the  local  velocity,  v,  and  local  depth,  d. 

e . Local  Froude  numbers  greater  than  3.0  should  be  avoided^; 
the  Froude  number  should  be  gradually  increasing  through- 
out the  transition. 

f^.  The  maximum  permissible  degree  of  curvature  of  the  side- 
walls  for  either  convergent  or  divergent  conditions  isA 


|Az/As|  = tan  *(— (1) 
3F 

— v — — 

where  F = — , — ; v and  d are  the  averages  of 

Vgd 

velocities  and  depths  at  the  beginning  and  at  the  end  of 
the  transition. 

A4 

£.  A convergence  should  be  such  that 


LAs  _ 1 

EAz  “ 0.382  - 0.116F 


(2) 


in  which  EAs  is  the  center-line  station  distance  mea- 
sured along  the  slope  downstream  from  the  center  line  of 
the  crest,  EAz  is  the  horizontal  convergence  of  one 
sidewall,  and  F is  the  local  Froude  number. 

9.  Surface  waves  in  supercritical  flow  are  analogous  to  shock 

FI  F2 

waves  in  compressible  gases.  ’ This  analogy  has  permitted  a reci- 

procity of  information  between  studies  concerned  with  liquid  flows  and 

those  with  gas  flows.  Two  early  works  that  demonstrate  this  information 

, F3 

exchange  are:  a study  by  Preiswerk,  who  applied  the  methods  of  gas 

F4 

dynamics  to  free  surface  water  flows;  and  a study  by  Laitone,  who  used 
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liquid  flows  in  a study  of  transonic  gas  dynamics.  A summary  of  the 

F4 

equations  which  describe  the  analogy  as  outlined  by  Crossley  is  as 
follows : 


Liquid  Flow  Gas  Dynamics 


in  which:  h = water  depth 

u , u = velocities  in  the  x and  y directions, 
xv  J 

J respectively 

g = acceleration  due  to  gravity 

T = absolute  temperature 

Y = ratio  of  specific  heats 

Cp  = specific  heat  at  constant  pressure 

The  analogous  quantities,  given  an  ideal  gas  having  y = 2,  are  gh  to 

2 2 

C T,  h/h  to  T/T  , h/h  to  p/p  , and  h /h  to  p/p  . Of  course,  the 
p o o o o o o 

analogy  is  subject  to  numerous  and  serious  limitations  (for  example, 
the  significance  of  energy  dissipation  as  well  as  the  physical  charac- 
teristics of  the  gas,  yt2,  must  be  taken  into  account).  For  more  de- 
tailed information  on  the  validity  of  this  analogy,  the  reader  is  re- 

F2  FI 

ferred  to  Crossley  and  to  Harleman. 

10.  The  analysis  which  is  the  current  basis  for  the  design  of 
transitions  for  supercritical  flow  is  summarized  in  a series  of  publi- 
cations presented  to  the  ASCE  Symposium  on  High-Velocity  Flow  in  Open 
Channels  in  1948.  In  contrast  to  the  empirical  observations  relevant  to 
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specific  design  problems  (some  are  listed  in  paragraph  8),  these  studies 
capitalized  on  the  insight  obtained  through  the  supercritical-supersonic 
analogy  described  above  to  obtain  the  basic  physical  principles  for  de- 
signing transitions  for  supercritical  flow. 

El 

11.  Ippen  describes  the  mechanics  of  supercritical  flow  for 
converging  boundaries.  The  practical  application  pertains  to  open- 
channel  flow  over  a horizontal  floor  in  which  surface  disturbances 
appear  as  a consequence  of  the  geometry  of  the  lateral  boundaries.  Two 
distinct  methods  of  approach  are  presented: 

a.  Gradual  surface  changes  are  analyzed  on  the  basis  of 
constant  specific  head. 

b.  Standing  wave  fronts  of  appreciable  height  (oblique  or 
slanting  hydraulic  jumps)  are  analyzed  with  a consider- 
ation of  the  energy  dissipation  involved. 

The  following  shock-wave  intersections  and  reflections  are  included: 

a.  Reflection  of  shock  wave  at  opposite  wall. 

jb.  Intersection  of  shock  waves  of  equal  intensity  and  of 
unequal  intensity. 

£.  Convergence  of  two  shock  fronts. 

d.  Shock  front  and  negative  deflection  waves. 

E2 

12.  Ippen  and  Dawson'  apply  the  mechanics  of  supercritical  flow 
directly  to  the  design  of  channel  contractions.  Examples  of  the  compu- 
tational procedure  are  presented,  and  comparisons  with  experimental  re- 
sults are  made.  Of  particular  interest,  experimental  water-surface 
elevation  contours  for  subcritical  (F<1)  and  supercritical  (F>1)  flows 
through  two  contractions  are  included. 

13.  High-velocity  flow  in  open-channel  expansions  is  analyzed  by 

F.3 

Rouse,  Bhoota,  and  Hsu  ; in  general,  water-surface  contours  were  deter- 
mined analytically  (by  means  of  the  method  of  characteristics)  for  nu- 
merous experimental  transition  designs  which  were  then  tested  in  the 
laboratory.  The  experimental  verification  is  reasonable  at  mild  slopes, 
less  so  at  steep  slopes  (10  percent  is  the  steepest  slope  tested).  For 
mild  slopes,  the  data  are  generalized  (not  a rigorous  identity)  in  the 
form  of 
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in  which  s and  z are  the  longitudinal  and  transverse  coordinates, 
respectively;  B is  the  channel  width;  and  the  subscript  1 refers 
to  conditions  at  the  entrance  to  the  expansion. 

14.  The  basic  formulation  and  design  procedures,  as  presented  by 
El 

Ippen  have  recently  been  incorporated  into  a computer  spillway  design 

C2 

procedure  by  Taubert.  The  computer  program  is  suited  to  funnel-shaped 
(straight  crest  in  Figure  1)  and  to  nozzle-shaped  (which  has  a reverse 
curve  in  the  sidewalls)  contractions;  the  program  is  not  applicable  to 
fan-shaped  contractions  (curved  crest  in  Figure  1) . The  program  con- 
sists of  five  basic  algorithms  as  follows: 

a.  Algorithm  1 is  used  to  calculate  the  characteristics  of 
standing  wave  fronts  caused  by  small  deflections  (specif- 
ic head  is  constant). 

b.  Algorithm  2 is  used  to  calculate  the  characteristics  of 
standing  waves  with  energy  dissipation  (oblique  shocks). 

£.  Algorithm  3 is  concerned  with  wave  reflection  from  a 

vertical  wall.  The  angle  of  deflection  at  the  reflected 
wave  front  is  set  equal  to  that  of  the  arriving  front; 
algorithm  1 or  2 (depending  on  energy  dissipation  con- 
siderations) is  used  to  advance  across  the  new  wave. 

cl.  Algorithm  4 is  concerned  with  the  intersection  of  two 
waves,  either  of  equal  or  of  unequal  intensity. 

e.  Algorithm  5 is  concerned  with  the  convergence  of  two  wave 
fronts. 

C2 

15.  Taubert  points  out  that  the  basic  equations  assume  a hori- 
zontal floor  and  uniform  flow  and  that  these  assumptions  are  not  valid 
in  practice.  These  two  effects  tend  to  compensate  for  one  another. 

However,  when  the  influence  of  the  floor  slope  is  strongly  dominant,  an 

E3 

elongation  of  the  flow  event  occurs.  The  effect  of  vertical  accelera- 
tion at  the  crest  of  a chute,  which  also  was  not  included  in  the  basic 

C3 

formulation,  is  described  by  Gunzel. 

16.  Most  flow  phenomena  occurring  in  steep  chutes  with  straight 
parallel  walls  are  also  of  concern  in  the  design  of  converging  chutes. 
Some  observations  and  studies  that  are  of  particular  concern  are  as 
follows: 


11 


a.  The  flow  near  the  crest  is  highly  nonunifonn;  however, 
the  upper  nappe  profile  in  this  region  may  be  approxi- 
mated by  means  of  successive  potential  flow  (flow  net) 
solutions.  The  uppermost  streamline  is  adjusted  between 
solutions  and  is  finally  located  such  that  the  total  head 
is  constant  along  the  free  surface  (see  HendersonA8  and 
Cassidy ,C4 ,C5  for  example).  Gunzel^  presents  experi- 
mental data  concerning  the  effects  of  the  rapidly  varied 
flow  near  the  crest  of  spillways.  Cox®l  presents  empir- 
ical curves  which  can  be  used  to  locate  the  upper  nappe 
profile  for  flows  over  crests  of  particular  shapes.  The 
slope  of  the  free  surface  at  control  points  is  the  sub- 
ject of  a recent  study  by  Wilkinson. ®2 

b.  The  gradually  varied  profile  farther  downstream  from  the 
crest  is  normally  an  S2  curve  which  may  be  calculated  by 
standard  methods  (see  Chow,A6  for  example) . The  USBR^-* 
design  guide  contains  a nomograph  that  simplifies  these 
backwater  computations  for  a steep  chute. 

c_.  Liggett  and  Vasudev®3  discuss  slope  and  frictional 

effects  in  two-dimensional  supercritical  flow.  A Corps 
of  Engineers  publication®^1  describes  the  development  of 
the  turbulent  boundary  layer  (and  a means  for  evaluating 
the  corresponding  energy  losses)  downstream  from  a spill- 
way crest. 

ci.  Air  entrainment  must  be  considered  for  flows  in  long 

chutes.  This  phenomenon  often  occurs  in  prototype  struc- 
tures; however,  it  is  a rare  occurrence  in  hydraulic 
models.  A recent  study  by  Keller,  et  al.,®5  discusses 
the  air  entrainment  phenomenon  and  presents  a means  for 
determining  the  point  along  the  chute  at  which  self- 
aeration of  the  flow  begins. 

e.  Flow  instability  can  occur  in  a steep  chute;  the  re- 
sulting unsteady  flow  is  termed  roll  waves  or  slug  flow. 
Roll  waves  are  discussed  in  numerous  publications^ >^2 ,G3 
the  presentation  by  MayerG3  and  the  discussions  thereof 
constitute  a reasonably  complete  background  of  the 
subject . 


Current  Study 

17.  Each  aspect  of  the  overall  spillway  design  problem  noted  pre- 
viously in  the  brief  literature  review  has  been  studied,  and  in  many 
regards  resolved,  by  earlier  investigations.  The  study  herein  is 


12 


therefore  reduced  to  the  following  items. 

a.  Construction  and  operation  of  a physical-model  test 
facility  to  provide  information  regarding: 

(1)  The  point  of  origin  along  the  wing  wall  at  which  the 
first  standing  wave  occurs  (convergence). 

(2)  The  geometry  of  the  first  wave  front  (convergence). 

(3)  The  geometry  of  the  flow  for  the  subsequent 
diverging-wall  situation  (divergence). 

b.  Suggesting  steps  to  be  taken  in  order  to  design  a 
satisfactory  shape  in  plan  of  a chute  sidewall. 


PART  II:  EXPERIMENTAL  APPARATUS 


Description 


18.  Figure  2 shows  the  experimental  facility;  details  of  the 
apparatus  are  presented  in  Plate  1. 


Figure  2.  Experimental  facility 
A summary  of  the  basic  components  follows: 

a.  The  supply  system,  from  the  sump  at  the  toe  of  the  chute 
to  the  headbay,  contained  a pumping  unit  made  up  of  four 
1-hp  centrifugal  pumps  in  parallel  that  discharged 
through  a 6-in.-diam  delivery  line.  The  maximum  capacity 
of  the  system  was  about  2.25  cfs;  flow  control  was  by 
means  of  a 6-in.  gate  valve  in  the  delivery  line. 

b.  The  dimensions  of  the  headbay  and  approach  channel  are 
shown  in  Plate  1.  The  geometry  is  such  that  the  flow 
lines  approaching  the  crest  are  parallel  ir.  plan.  Horse- 
hair and  screen  matting  were  used  as  baffles  in  the  head- 
bay  so  that  surface  irregularities  were  essentially 
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eliminated  before  reaching  the  approach  channel. 

£.  The  dimensions  of  the  chute  spillway  are  also  shown  in 
Plate  1.  Although  the  slope  of  the  chute  was  variable, 
a specific  slope  of  0.1799  on  1 (equaling  10.2  degrees) 
was  used  in  all  tests  reported  herein.  Details  of  the 
profile  of  the  spillway  crest  are  shown  in  Figure  3. 


Figure  3.  Details  of  the  spillway  crest 

d.  The  basic  measuring  equipment  consisted  of  a 90°  bend 
meter  for  discharge  and  point  gages  for  water-surface 
elevation.  Dye  and  confetti  were  used  on  occasion  for 
visual  observations  of  flow  patterns. 


General  Hydraulic  Parameters 


19.  The  three  basic  coefficients  that  typify  the  apparatus  in- 
volved (a)  calibration  of  the  bend  meter,  (b)  determination  of  a free- 
flow crest  discharge  coefficient,  and  (c)  evaluation  of  the  hydraulic 
roughness  of  the  steep  chute.  These  items  are  discussed  in  the  follow- 
ing paragraphs. 

20.  An  approximate  relationship  for  discharge  measurement  by 
means  of  a 90°  bend  meter  is 


a/2gAh 


(4) 
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in  which  r is  the  center-line  radius  of  the  pipe  bend,  D and  a are 
respectively  the  pipe  diameter  and  cross-sectional  area,  g is  the 
gravitational  acceleration,  and  h is  the  difference  in  piezometric 
head  between  the  inside  and  outside  of  the  bend.  For  the  experimental 
setup  shown  in  Plate  1,  r was  0.667  ft,  D was  0.5  ft,  and  a was 
0.1963  ft2. 

21.  The  calibration  equation  for  the  meter  was  established  in 
place  by  means  of  two  pitot-tube  traverses  across  the  outlet;  the  equa- 
tion of  the  meter  was  established  in  the  form 

Q - CM  (5) 

in  which  Cw  is  a coefficient  and  R is  the  manometer  differential  in 
inches  of  Meriam  M-3  fluid  (specific  gravity  equals  2.95  at  60°F) . 
Equation  4 yields  a value  of  Cu  for  use  in  Equation  5 equaling  about 
0.52;  the  in  place  average  calibration  value,  CM  equals  0.55,  was  used 
in  all  data  reduction  herein  and  should  be  accurate  within  ±5  percent 
over  the  range  of  discharges  of  interest.  Both  curves  are  shown  in 
Figure  4;  the  Reynolds*  number,  Rg  , ranged  from  1.4(10^)  to  4.9(10^) 
for  the  discharges  used  in  the  test  program. 

22.  The  crest  discharge  coefficient,  C , was  evaluated  from  the 
following  equation^ 

Q • C L'  H 3/2  (6) 

e 

in  which  L'  is  the  effective  crest  length  and  H , which  is  the  ele- 

e 

vation  of  the  energy  grade  line  of  the  approach  flow  relative  to  the 

crest,  equals  h + h as  previously  shown  in  Figure  1.  For  a crest 

V A1 

with  one  rectangular  abutment 

L’  - L - 0.1  H (7) 

c e 

— 

* Re  ■ — in  which  V - Q/a  is  the  mean  velocity  in  the  delivery  line 
and  v is  the  kinematic  viscosity. 
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RATE  OF  FLOW,  CFS 
Figure  4.  Bend  meter  calibration 

in  which  L is  the  crest  length  (2.500  ft  for  the  test  facility).  A 
c 

series  of  observations,  tests  H1-H14,  were  performed  in  order  to  deter- 
mine C as  a function  of  H : the  data  and  calculations  are  tabulated 

e 

below,  and  the  results  are  summarized  in  Figure  5. 


Test 

Q 

h 

h 

V 

H 

e 

L’ 

No. 

cfs 

ft 

ft 

ft 

ft 

C 

HI 

0.55 

0.143 

0.001 

0.144 

2.486 

4.05 

H2 

0.78 

0.181 

0.002 

0.183 

2.482 

4.01 

H3 

0.95 

0.207 

0.003 

0.210 

2.479 

3.98 

H4 

1.10 

0.229 

0.004 

0.233 

2.477 

3.95 

H5 

1.23 

0.248 

0.005 

0.253 

2.475 

3.91 

(Continued) 
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Test  Q h hv  He  L' 


No. 

cf  s 

ft 

ft 

ft 

ft 

C 

H6 

1.35 

0.262 

0.005 

0.267 

2.473 

3.96 

H7 

1.46 

0.274 

0.006 

0.280 

2.472 

3.99 

H8 

1.56 

0.289 

0.006 

0.295 

2.471 

3.94 

H9 

1.65 

0.299 

0.007 

0.306 

2.469 

3.95 

H10 

1.74 

0.309 

0.007 

0.316 

2.468 

3.97 

Hll 

1.82 

0.319 

0.008 

0.327 

2.467 

3.95 

H12 

1.91 

0.328 

0.008 

0.336 

2.466 

3.98 

HI  3 

1.98 

0.335 

0.009 

0.344 

2.466 

3.98 

H14 

2.01 

0.338 

0.009 

0.347 

2.465 

3.99 

c 


Figure  5.  Crest  discharge  coefficient 
(parallel  sidewalls) 

23.  The  average  value,  C = 3.97,  is  used  herein.  The  largest 
difference  between  an  observed  value  and  the  mean  is  2 percent,  which  is 
substantially  less  than  the  ±5  percent  confidence  in  the  discharge 
measurement . 

24.  A sketch  of  the  flow  profile  downstream  from  the  crest  is 
presented  in  Figure  6.  The  two  parameters  of  interest  are  the  ratio  of 
depth  at  the  crest  to  critical  depth,  ^O.OO^c’  anc*  t*le  evaluation  of 
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RAPIDLY  VARICD  FLOW 


Figure  6.  General  flow  profile  (parallel  sidewalls) 
the  normal  depth,  d^,  as  a function  of  unit  discharge,  q.  Data  re- 
garding the  depth  at  the  crest  are  tabulated  below  and  shown  in  Figure  7; 

A9 

the  critical  depth  is  evaluated  using  the  equation 


d = 

c 


(8) 


Test 

No. 

q 

cfs/f t 

d 

c 

ft 

do.oo 

ft 

do.oo 

d 

c 

HI 

0.221 

0.115 

0.108 

0.939 

H2 

0.314 

0.145 

0.139 

0.958 

H3 

0.383 

0.166 

0.163 

0.978 

H4 

0.444 

0.183 

0.179 

0.979 

H5 

0.497 

0.197 

0.194 

0.985 

H6 

0.546 

0.210 

0.209 

0.996 

H7 

0.591 

0.221 

0.219 

0.991 

H8 

0.631 

0.231 

0.232 

1.001 

H9 

0.668 

0.240 

0.242 

1.006 

H10 

0.  705 

0.249 

0.248 

0.994 

Hll 

0.738 

0.257 

0.256 

0.997 

H12 

0.775 

0.265 

0.275 

1.000 

HI  3 

0.803 

0.272 

0.271 

0.997 

H14 

0.815 

0.274 

0.277 

1.008 

25.  Five  measured  water-surface  profiles  (tests  NH1-NH5)  are 
shown  in  Figure  8.  Because  of  the  limited  length  of  the  chute,  the 
flow  became  essentially  uniform  only  for  small  discharges.  For  example. 
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Figure  7.  Depth  at  the  spillway  crest 
(parallel  sidewalls) 


in  Figure  8,  only  the  lowest  flow  rate,  0.319  cfs/ft,  appears  to  be  near 
a uniform  flow  condition  at  the  extreme  measurement  location,  sta 
2+0.50  ft.  This  test  condition  was  used  to  evaluate  a roughness  coef- 
ficient, n,  for  Manning's  flow  equation 


n 


1.486 


q 


(9) 


With  normal  depth  d * 0.049  ft,  S.  = S * 0.1771,  and  q = 0.319 

n f o 

cf 8/ ft  in  Equation  9,  a value  of  n * 0.0124  was  obtained.  The  value 

A6 

n « 0.012,  which  agrees  with  published  values  for  a planed  wood  sur- 
face, was  used  in  subsequent  computations  herein. 

26.  S2  backwater  profiles  for  the  five  flows  are  also  shown  in 
Figure  8;  these  profiles  are  computed  numerically  in  the  form 

y1+1  - yt  + Ay  (10a) 


20 


Figure  8.  S2  backwater  profiles  (parallel  sidewalls) 


i+1 


-1  (V 

x +2g_^i±l_ 
i S 


2 2n 
Vi  > 


+ yi+l  ~yi 


o - Sf 


(10b) 


A definition  sketch  is  shown  in  Figure  9.  In  performing  the  numerical 

integrations  y is  held  constant,  y^  “ yc , x^  = 0.0,  and  S^.  is 

evaluated  using  Manning's  equation.  In  plotting  Figure  8,  the  slope 

coordinate,  s,  was  used  rather  than  the  horizontal  coordinate,  x, 

and  the  backwater  profiles  are  arbitrarily  offset  2.00  ft  along  the  s- 

axis.  Because  of  the  rapidly  varying  flow  conditions  near  the  crest  the 

S2-profiles  and  the  measured  profiles  are  dissimilar  except  at  larger 

A5 

s values;  consequently,  the  computational  technique,  wherein  the  off- 
set at  the  crest  is  used  to  obtain  closer  agreement  between  calculated 
and  actual  water-surface  profiles,  is  not  well  suited  to  the  profile 
near  the  crest. 


a 
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f 


Figure  9.  Definition  sketch  (gradually  varied  flow) 

27.  The  amount,  L , by  which  the  straight  wing-wall  segments 
w 

must  be  extended  upstream  in  order  to  obtain  parallel  flow  at  the  crest 
was  determined  by  a separate  series  of  tests.  Eight  wing-wall  segments, 
ranging  from  0.33  to  3.0  ft,  were  used.  The  discharge  was  gradually 
increased  until  the  nonparallel  streamlines  (as  observed  by  means  of  dye 
streaks  in  the  flow)  originating  from  the  vortex  at  the  leading  edge  of 
the  walls  persisted  to  the  crest.  The  two  photographs  in  Figure  10  show 
the  vortex  with  two  wing-wall  conditions  at  a low  discharge  value.  The 
results  indicated  that  the  minimum  acceptable  value  for  the  complete 
range  of  discharges  studied  in  the  test  chute  was  1.33  ft  or  iWp 
equals  3.8;  this  extension  was  used  for  all  subsequent  tests. 
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PART  III:  TEST  OBSERVATIONS  AND  RESULTS 


28.  The  following  observations  are  presented  on  the  basis  of  pre- 
liminary visual  inspection  of  various  flows  in  the  chute  coupled  with 
the  previous  literature  survey. 

a.  Approach  conditions.  The  comments  noted  previously 
(paragraph  3)  regarding  the  need  for  a uniform  approach 
flow  aligned  normal  to  the  crest  are  again  emphasized. 

Any  obstruction  or  misalignment  not  subject  to  prior 
study  may  result  in  unpredicted  undesirable  waves  in  the 
chute;  suspicion  of  this  situation  warrants  a hydraulic 
model  study  during  the  design  phase.  The  hydraulic  model 
currently  is  the  only  proven  means  of  modifying  a design 
in  order  to  cope  with  poor  approach  flows. 

b.  Wing -wall  extensions.  The  current  study  used  straight 
wing  walls;  these  walls  must  extend  a substantial  dis- 
tance (3.8p  minimum  for  the  conditions  studied)  upstream 
from  the  crest  in  order  to  reduce  the  effects  of  the  up- 
stream vortices.  In  practice,  these  walls  are  usually 
curved  in  plan,  thereby  reducing  the  upstream  vortices 
and  making  shorter  wing  walls  feasible. 

£.  Rapidly  varied  flow  at  the  crest.  The  gradually  varied 
S2  flow  profile  does  not  approximate  the  flow  conditions 
immediately  downstream  from  the  crest  even  when  parallel 
sidewalls  are  used  (see  Figure  9).  Hence,  either  hydrau- 
lic model  data  or  a numerical  solution  that  takes  into 
account  the  rapid  acceleration  at  the  crest  is  needed. 
Model  data  are  available  for  most  standardized  crest 
shapes . 

29.  The  following  observations  have  to  do  with  the  effects  on 
the  flow  due  to  the  convergence  and  divergence  of  the  chute  sidewalls 
immediately  downstream  from  the  crest.  Plates  2,  3,  and  4 show  the 
cross-sectional  flow  profiles  for  a straight  wall  converging  section; 
Plate  5 shows  a straight  wall  converging  followed  by  a straight  wall  di- 
verging section.  The  s and  z coordinate  system  is  also  shown  in 
Plate  5.  The  data  indicate  the  following. 

£.  Flow  profile  (cross  section)  at  the  crest.  With  straight 
approach  flow,  a two-dimensional  flow  situation  exists 
and  the  depth  at  the  crest  is  approximately  equal  to 
critical  depth  as  previously  shown  in  Figure  8.  On  the 
other  hand,  with  a converging  wing  wall,  the  flow  is 
three-dimensional,  for  which  the  depth  at  the  crest  in- 
creases nearer  to  the  wall.  Evaluating  critical  depth. 
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the  observed  change  is 


dc  , as  dc  = 
as  follows: 


Q 

p/H 

e 

d0.000/dc 

Test 

cfs 

Near  Wing  Wall 

z = 0.0 

Nl 

0.590 

2.29 

1.13 

0.99 

N2 

0.992 

1.62 

1.09 

0.98 

N3 

1.481 

1.22 

1.05 

0.96 

! Of 

standing  wave  (convergence) . If  the 

design  head, 

for 

the  crest 

is  assumed 

to  be  twice  the 

radius  of 

tRe  crest,  then  the  values  of  H /H  for  tests  Nl,  N2, 
and  N3  are  0.73,  1.032  and  1.373^  respectively.  The  data 
in  Plates  2-4  indicate  that  for  these  flow  conditions  the 
initial  standing  wave  has  a "gradual"  surface  change 
rather  than  being  a standing  wave  of  "appreciable"  height 
even  with  a maximum  recommended  wing-wall  convergence 
(|As/Az|  =5).  In  addition,  even  at  the  maximum  attain- 
able flow  rates,  H /H„  = 1.67,  no  well-defined  or  sud- 
e D 

den  change  in  water  depth  occurred. 

£.  Spreading  of  the  initial  standing  wave  across  the  chute. 

A small  gravity  wave  tends  to  move  across  the  chute 
according  toEi 

6 * sin_1(j)  (11) 

in  which  E is  the  angle  of  the  wave  and  F = v/*/gd  is 
the  local  Froude  number.  A definition  sketch  and  a plot 
of  Equation  11  are  given  in  Plate  6.  The  calculated 
variation  of  6,  using  the  recorded  depths  along 
z = 0.0  in  Equation  11,  are  shown  in  Figure  11.  Data 
indicate  that  in  the  vicinity  of  the  crest,  the  wave 
angle  can  be  represented  by  a single  function, 

3 = / (s/H  ),  for  the  test  chute.  Using  this  function 
and  integrating  numerically  along  the  chute 


1/3 


[(3.97)2/g]  H 


s 


i+1 


z 


i+1 


s^  + As 

A + ^ 

z^  + As  tan  ( 2 ) 


(12a) 

(12b) 


gives  the  spread  of  the  wave  for  tests  Nl,  N2 , and  N3 
shown  in  Figure  12.  The  initial  conditions  (at  i = 0) 
are : 


s = 0.5H 
o e 

6 = 85° 

o 


(13a) 

(13b) 
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(13c) 


Figure  11.  Orientation  of  a small  gravity  wave  (z  = 0.0) 

z - 0.940  - s tan  10.2° 
o o 


Locations  of  the  computed  wave  are  also  shown  in  Plates 
2,  3,  and  4 and  are  in  general  agreement  with  the  mea- 
sured profiles. 


d. 


Depths  along  the  wall.  For  gradual  wall  deflections  into 
the  flow^l 


0+0 


1 


/r 


tan 


V 


f2-i 


(14) 


in  which  0 is  an  initial  condition  (F  * F^) , and 
0+0.  is  the  angle  of  deflection.  Equation  14,  with 
1 <_  F~<  4 is  also  plotted  in  Plate  6;  the  value  of 
0+01  varies  from  0°  for  an  infinite  Froude  number  to 
65.88°  for  F = 1.  Data  corresponding  to  Equation  14  are 


Z,  FT 


Figure  12.  Spreading  of  the  initial  standing  wave 

plotted  in  Figure  13.  Using  the  values  of  F along 
z = 0.0,  obtained  from  Figure  11,  values  of  (up- 
stream from  the  disturbed  flow)  are  calculated.  Since 
the  overall  deflection  0 = 11.31°  is  known,  the  values 
of  0+0^  are  also  plotted.  The  experimental  data 
(0  + 0^)  are  evaluated  using  measured  depths  along  the 
wall  to  evaluate  F2  for  tests  Nl,  N2,  and  N3  and  are 
also  shown  in  Figure  13.  These  data  are  in  reasonable 
agreement  with  the  values  derived  directly  from  Equation 
14  using  the  mainstream  F^. 

e.  Surface  contours  in  diverging  flow.  Channel  expansions 
are  analyzed  by  Ippen^^J in  general,  in  expanding  flow 
no  steep  wave  fronts  exist  and  the  two  following 


27 


Figure  13.  Conditions  along  a converging  sidewall  (0  = 11.31°) 
conditions  must  be  satisfied: 

(1)  There  is  an  essentially  constant  specific  head  for 
the  geometric  extent  of  the  flow  considered. 

(2)  The  wave  angle,  B,  is  equal  everywhere  to  the  local 
value  of  sin-*  (1/F). 

The  surface  contours  of  the  entire  field  of  flow  are  then 
mapped  via  the  method  of  characteristics.  Rouse,  Bhoota, 
and  Hsue^  conducted  experiments  on  channel  expansions  and 
found  that  the  surface  contours  were  adequately  described 
by  a simplified  empirical  function  (Equation  3 in  para- 
graph 13).  Their  experimentally  determined  surface  con- 
tours for  a 90°  expansion  on  a horizontal  floor  with  Fj 


28 


0 0.5  1.0  1.8  2.0 


S-Si 
B,  F, 

Figure  14.  Divergence  (surface  contours) 

equals  2.0  are  approximately  as  shown  in  Figure  14.  Five 
comparable  data  values  for  the  11.31°  expansion  (Tests  N2 
and  N3  are  also  shown  in  Figure  14.  For  these  tests,  the 
value  of  Bi  is  assumed  equal  to  the  width  of  the  stand- 
ing wave  at  the  end  of  the  convergence  as  obtained  from 
Figure  12.  In  general.  Rouse,  Bhoota,  and  Hsu^  have 
previously  established  a reasonably  complete  data  set  for 
expansions  in  supercritical  flow  and  since  they  have  ver- 
ified the  accuracy  of  the  method  of  characteristics  solu- 
tion for  expansions,  this  aspect  of  the  problem  is  not 
pursued  further  herein.  The  current  data  tend  to  indi- 
cate that  the  previous  work  can  be  satisfactorily  applied 
to  supercritical  flow  in  a steep  chute  having  convex 
sidewalls  provided  the  divergence  is  not  in  a highly  non- 
uniform  flow  region. 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


30.  The  brief  literature  survey  represents  a very  small  part  of 
the  published  literature  relevant  to  flow  in  steep  chutes.  However,  the 
survey,  coupled  with  the  limited  experimental  data  herein,  indicate  the 
following: 

a.  The  author  can  conceive  of  no  sidewall  configuration 
(other  than  parallel  walls)  that  will  result  in  a hori- 
zontal water  surface  across  the  chute.  For  this  reason, 
it  is  suggested  that  the  designer  evaluate  the  water- 
surface  contours  for  each  proposed  design  and  reject 
those  plans  that  result  in  unacceptable  water-surface 
variations . 

b.  Given  a satisfactory  approach  condition  (i.e.,  uniform 
flow  normal  to  the  crest  and  properly  designed  wing 
walls)  and  using  state-of-the-art  information  on  chute 
design,  a computer  model  for  spillway  design  can  be  de- 
veloped. This  model  would  take  into  account  the  varia- 
tion of  the  flow  cross  section  along  the  crest  (para- 
graph 29a),  the  initial  convergence  (paragraph  29b),  the 
effect  of  rapidly  varied  flow  on  the  spreading  of  the 
standing  wave  (paragraph  29c),  and  the  increased  depths 
along  the  wall  in  the  converging  section  (paragraph  29d). 
The  surface  contours  for  tne  flow  along  the  subsequent 
diverging  wall  (paragraph  29e)  can  be  established  using  a 
programmed  method-of-characterist ics  solution.  A model 
that  neglects  any  of  the  above  factors  cannot  provide  a 
solution  that  is  entirely  satisfactory  to  the  designer. 

£.  Given  an  unsatisfactory  approach  condition  it  is  unlikely 
that  a computer  model  can  provide  a reliable  sidewall 
design.  For  this  particular  circumstance,  the  hydraulic 
model  is  probably  the  best  available  design  tool. 
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APPENDIX  A:  NOTATION 


DIMENSION 


SYMBOL  DESCRIPTION (F-L-T) 


a 

cross-sectional  area 

L2 

B 

approach  width  (supercritical  flow) 

L 

B 

approach  width  (subcritical  flow) 

L 

C 

crest  discharge  coefficient 

il'V1 

CM 

bend  meter  coefficient 

l5/2t-‘ 

C 

specific  heat  at  constant  pressure 

- 

P 

(Btu/lb  °R) 
m 

d 

local  depth  (normal  to  bottom) 

L 

D 

diameter 

L 

d 

c 

critical  depth 

L 

d 

n 

normal  depth 

L 

1 

i 

d0.00,d0.71’ 

depth  at  sta  0.00,  0.71,  etc. 

L 

f 

F 

Froude  number  (F  = v/i/gd) 

none 

g 

gravitational  acceleration 

LT-2 

h 

piezometric  head 

L 

h 

V 

velocity  head 

L 

H 

e 

head  on  crest  (H  = h + h ) 
e v 

L 

hd 

design  head 

L 

i 

integer  subscript 

none 

L 

c 

crest  length 

L 

L' 

effective  crest  length 

L 

Lw 

wing-wall  extension  length 

L 

n 

Manning's  roughness  coefficient 

l"1/3t 

p 

crest  height 

L 

q 

unit  discharge 

lV1 

Q 

total  discharge 

lV1 

r 

radius 

L 

R 

manometer  differential 

L 

R 

e 

VD 

Reynolds  number  (R^  * — ) 

none 

s 

dimension  (measured  parallel  to 

L 

channel  bottom) 

A1 

v 

— - 

SYMBOL 

DESCRIPTION 

DIMENSION 

(F-L-T) 

S 

slope 

none 

T 

temperature  (°R) 

- 

C 

X 

c 

velocity  components 

V 

local  velocity 

LT 

V 

mean  velocity 

LT-1 

X 

dimension  (horizontal) 

L 

y 

dimension  (vertical) 

L 

z 

dimension  (transverse) 

L 

B 

deflection  angle  (wave) 

none 

e 

deflection  angle  (wall) 

none 

Y 

ratio  of  specific  heats 

none 

V 

kinematic  viscosity 

lV1 

P 

density 

fl'4t2 

Sf 

slope  of  the  energy  grade  line 

none 

s 

o 

slope  of  the  channel  bottom 

none 

A2 
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